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We analyze the translational and rotational motion of an ellipsoidal Brownian particle from the 
viewpoint of stochastic thermodynamics. The particle’s Brownian motion is driven by external 
forces and torques and takes place in an heterogeneous thermal environment where friction coeffi¬ 
cients and (local) temperature depend on space and time. Our analysis of the particle’s stochastic 
thermodynamics is based on the entropy production associated with single particle trajectories. 
It is motivated by the recent discovery that the overdamped limit of vanishing inertia effects (as 
compared to viscous fricion) produces a so-called “anomalous” contribution to the entropy produc¬ 
tion, which has no counterpart in the overdamped approximation, when inertia effects are simply 
discarded. Here, we show that rotational Brownian motion in the overdamped limit generates an 
additional contribution to the “anomalous” entropy. We calculate its specific form by performing 
a systematic singular perturbation analysis for the generating function of the entropy production. 
As a side result, we also obtain the (well-known) equations of motion in the overdamped limit. We 
furthermore investigate the effects of particle shape and give explicit expressions of the “anomalous 
entropy” for prolate and oblate spheroids and for near-spherical Brownian particles. 


I. INTRODUCTION 

The theory of Brownian motion m, developed in dif¬ 
ferent formulations by Einstein [4] , Smoluchowski [5] and 
Langevin [5] around 1905 and 1906, describes the dynam¬ 
ics of a particle suspended in a fluid. A prototypical ex¬ 
ample is a small colloidal object, e.g. a polystyrene bead 
about a micrometer in size, floating in water at room tem¬ 
perature. Even without the action of externally applied 
forces, the particle is in an animated and erratic state 
of motion, generated on microscopic scales by collisions 
with the water molecules and visible on mesoscopic scales 
as irregular diffusive movement. Based upon the works 
mentioned above, this Brownian motion is most success¬ 
fully modelled (on the mesoscopic scales) by stochastic 
differential equations CHa], augmenting the Newtonian 
equations of motion for the particle by the forces from 
the surrounding liquid, namely Stokesian friction propor¬ 
tional to the particle velocity and thermal fluctuations re¬ 
lated to the fluid temperature. Both originating from the 
surrounding fluid bath as their source, the strengths of 
these two forces are related by the fluctuation-dissipation 
theorem laiiii. 

Typically, this set of equations of motion for parti¬ 
cle position and velocity can be simplified by adopting 
the so-called overdamped approximation, where one com¬ 
pletely neglects inertia effects in the particle dynamics. 
This procedure is justified because for micrometer-sized 
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objects suspended in water friction forces are by orders 
of magnitude larger than inertial forces [iniiii], so that 
on time- and length-scales accessible under typical exper¬ 
imental conditions only the overdamped behavior is ob¬ 
servable. Indeed, it takes a delicate experimental effort 
to actually measure the velocity of a micrometer-sized 
Brownian particle in a liquid [mils]. 

In recent years it has been demonstrated that the 
stochastic equations of motion not only describe the ir¬ 
regular particle dynamics to high accuracy, but are also 
a valid starting point for a consistent theory of thermo¬ 
dynamic quantities which are associated with the parti¬ 
cle movement, such as heat, work or entropy production 
[HHni, even if the particle is driven away from ther¬ 
mal equilibrium conditions with the heat bath. In this 
newly emerging field called stochastic thermodynamics, 
central results of surprising generality have been discov¬ 
ered in form of fluctuation relations; a recent summary is 
provided in m- The entropy production (rate) plays a 
particularly important role, because it constitutes a mea¬ 
sure of irreversibility by relating probabilities of particle 
trajectories to their time-reversed counterparts |18j . 

In the present context of a colloidal particle whose 
stochastic dynamics is well-described by overdamped 
equations of motion, the fact that the definition of such a 
central concept as entropy production (rate) is based on 
single trajectories immediately raises the question of how 
neglecting the velocity degrees of freedom may affect the 
particle’s stochastic thermodynamics. Investigating this 
question, it has been shown in [19] that the overdamped 
approximation does not fully capture the entropy produc¬ 
tion rate if the thermal environment is inhomogeneous 
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(but in equilibrium locally). Rather, there is a contribu¬ 
tion to the entropy production in addition to those pre¬ 
dicted from the overdamped approximation, which can 
not be obtained from the statistics of the overdamped 
trajectories. Being connected to the breaking of time and 
velocity reversal symmetry, this phenomenon has been 
dubbed “entropic anomaly” |19] , in analogy with similar 
anomalies encountered in physics. 

Perhaps the best known such anomaly is the viscous 
dissipative anomaly in turbulence: the energy dissipation 
rate in the fluid remains finite in the limit of vanishing 
viscosity, while it vanishes for the viscosity being set to 
exactly zero [2^. This anomaly is caused by the break¬ 
ing of time-reversal symmetry; time-reversal symmetry is 
present under strictly inviscid conditions, but not for fi¬ 
nite, yet arbitrarily small viscosity. Similarly, in quantum 
physics an anomaly is related to a symmetry operative on 
the classical level (Planck’s constant set to zero), which 
is broken in the corresponding quantum theory mm- 

The occurrence of the “anomalous entropy production” 
in the overdamped limit has been discovered in [19] for 
purely translational motion of a spherical Brownian par¬ 
ticle through an inhomogeneous thermal environment, 
see also [53]. It has been further analyzed in |21| and 
shown to occur in general classes of stochastic systems, 
including discrete stochastic processes; a related analy¬ 
sis of discrete processes is presented in [25 . The influ¬ 
ence of a gradient flow on the “entropic anomaly” has 
been studied in |26|, where also effects of particle rota¬ 
tion have been taken into account. Moreover, the “en¬ 
tropic anomaly” has been demonstrated to affect opti¬ 
mal stochastic transport, i.e. driven processes which op¬ 
timize entropy production m. to induce an efficiency 
loss in microscopic stochastic heat-engines [2H] , and even 
to play a non-trivial role in microevolution |29j . Finally, 
the “entropic anomaly” can be seen as an explicit exam¬ 
ple of the general observation that the entropy produc¬ 
tion may depend on the scale of description [2il[3nil3Il; 
i.e. coarse-graining the dynamical equations of a physical 
system by integrating out a sub-set of degrees of freedom 
typically reduces entropy production 1301 m- Further 
concrete systems with such scale-dependent entropy pro¬ 
duction are a harmonic chain of two Brownian particles 
in contact with two different heat baths, having finite 
entropy production which is reduced to zero when one 
of the Brownian oscillators is integrated out [32] j and a 
dimer consisting of two Brownian particles at different 
temperatures with a harmonic (but stiff) vs. a rigid cou¬ 
pling [33] . 

In the present paper, we analyze in detail the Brownian 
motion of a non-spherical particle, especially the contri¬ 
butions of rotational degrees of freedom to the stochastic 
entropy production and to the appearance of “anomalous 
entropy” in the overdamped limit. We restrict ourselves 
to the case where there is no hydrodynamic coupling be¬ 
tween rotational and translational degrees of freedom. 
Our theory is thus valid for any particle with three mu¬ 
tually perpendicular symmetry planes |34j . including, in 


particular, the large class of spheroids and ellipsoids, but 
also rods and other rod-like shaped objects. 

This paper is organized as follows. In Sec. II, we 
describe the model, provide the fundamental governing 
equations for translational and rotational Brownian mo¬ 
tion, and introduce the basic path-wise thermodynamic 
quantities, such as heat and entropy production, follow¬ 
ing the approach of stochastic thermodynamics. Section 
III presents the derivation of the overdamped limit for the 
concepts introduced in Sec. II. Sections IV and V discuss 
the resulting overdamped dynamics and overdamped en¬ 
tropy production, including the anomalous contribution. 
In Sec. VI, the anomalous entropy production is explic¬ 
itly calculated for prolate and oblate spheroids. Sec. VII 
treats the case of slightly deformed spherical particles. 
We conclude with a short summary and discussion in 
Sec. VIII. The Appendices A-C contain some additional 
information or details of the calculations. 


II. DYNAMICS AND ENTROPY PRODUCTION 

The dynamics of the particle is governed by external 
forces and torques, thermal fluctuations, and viscous fric¬ 
tion. Our starting point is to take into account inertia 
effects as well, so that we describe such driven Brown¬ 
ian motion by a set of Langevin-Kramers equations [7H2] 
for the particle’s translational and rotational degrees of 
freedom. The translational motion of the center of mass 
X = (xi, X 2 , X 3 ) of the particle (with mass m) and its ve¬ 
locity V = (ui, W 2 , U 3 ) is modeled in the laboratory frame, 

X = V , (la) 

mi) = —71? -I- / -l- , (lb) 

where / = (/i,/ 2 ,/ 3 ) summarizes all deterministic ex¬ 
ternal forces, T is the temperature (/cb Boltzmann’s 
constant), and ^(t) = (Ci(t),^2(^)1 ^3(^)) unbiased 
Gaussian noise sources with correlations {^i{t)^j{t')) = 
6 ijS(t — t'). Finally, 7 is the translational viscous fric¬ 
tion tensor of the ellipsoid. This 3x3 tensor is positive 
definite and symmetric, so that it has a unique square 
root, which we represent by 7^^^, meaning (71/2)^^172 _ 
^i/2^i/2 _ ry symbol^ labels the matrix transpose). 

While the translational dynamics is represented in the 
laboratory frame, it turns out to be more convenient to 
write the rotational motion in a body-fixed frame (with 
origin in the particle’s center of mass), because then the 
inertia tensor / is independent of particle orientation and 
thus constant in time. The rotational Langevin-Kramers 
equation for the angular velocity u) = (wi,a; 2 ,wa) of 
the ellipsoid is then given by Euler’s equation of rigid 
body dynamics |35j with a total torque including not only 
externally applied torques but also viscous friction and 
thermal noise, 

lu) -I- a; X (/oi) = —rju) -|- M -t- y/2kBTr]^^'^C{t). (2a) 
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Here, ry is the 3x3 symmetric and positive definite ro¬ 
tational friction tensor with a unique square root 
M = {Ml, M2, M3) represents the deterministic torques 
acting on the particle, and C{t) = (Ci(0) C2(i)) C3(0) 
are unbiased and delta-correlated Gaussian noise source, 
which are independent of the translational ones ^{t). 

In order to specify the orientational position of the 
ellipsoid uniquely we choose two orthogonal unit vectors 
n = {ni,n2,ri3) and m = {mi, m2, m3), which are rigidly 
attached to the particle (see Fig. [^. Their movement is 
dictated by the angular velocity according to the kine¬ 
matic equations 



FIG. 1. Sketch of an ellipsoidal Brownian particle. The two 
vectors n and m are attached to the particle to parametrize 
its rotational position, see also main text, in particular 
Eqs. (2bl. 


n = LJ X n , rh = u} X m . (2b) 


Since the lengths of these two vectors are set to unity, 
and their relative orientation is kept fixed (we choose 
n-m — 0 ), this representation of the particle orientation 
has three free parameters, as expected for a represen¬ 
tation of rotation in three dimensions. The results we 
derive in this paper are of course independent of the spe¬ 
cific representation of the particle rotation. Instead of 
(2b) one could choose a quaternion representation [5S] . 
an Euler angle representation [53 [3S] , and also a differ¬ 
ential geometric representation in terms of local charts, 
as discussed in [53. For the calculation we are going 
to perform it turns out, however, that among the global 
representations the one in (2b) is the most convenient. 


We emphasize again that translational motion ([^ is 
written in the laboratory frame, while rotational dynam¬ 
ics is given in a reference frame fixed to the parti¬ 
cle with the origin being located at the center of mass 
of the particle. We do not distinguish quantities in the 
different reference frames by explicit labels though, in 
order to keep notation as simple as possible. In other 
words, throughout the present manuscript we follow the 
rule that all quantities associated with translation are 
represented in the laboratory system, while quantities as¬ 
sociated with rotation are represented in the body fixed 
coordinate system. 


A central property of our setup is the presence of an 
inhomogeneous thermal environment, i.e. we allow the 
temperature T to depend on space and time, but with 
thermal equilibrium being valid locally. Likewise, the 
friction tensors 7 and 77 are assumed to be functions of 
space and time, for instance because the fluid viscos¬ 
ity changes with the spatial variations of temperature 
or due to hydrodynamic effects close to boundaries m- 
Moreover, since translation is considered in the labora¬ 
tory frame, the translational friction tensor depends on 
the particle orientation, while the body-fixed rotational 
friction tensor does not. Concerning the external deter¬ 
ministic forces and torques, we allow for the most gen¬ 
eral case and take into account variations on position, 
orientation and time. Therefore, we have the following 
functional dependencies for the quantities appearing in 


( 0 ,( 121 ): 


T = T{x,t), 

(3a) 

7 = -f{x,n,m,t), 

(3b) 

V = V{x,t), 

(3c) 

f = f{x,n,m,t ), 

(3d) 

M = M{x, n, m, t). 

(3e) 


The equations of motion 0. 0 generate, for given 
initial position and orientation of the ellipsoid, stochas¬ 
tic trajectories in translational and orientational config¬ 
uration space. According to stochastic thermodynamics 
[m US], an entropy production AS is associated with 
such stochastic trajectories, which involves contributions 
from the change of particle entropy ASp and from the 
entropy production in the environment ASenv, 

AS = ASp + ASenv . (4) 

We briefly summarize the well-known results for transla¬ 
tional Brownian motion [38] and extend them to include 
rotational stochastic dynamics. 

The particle entropy is defined as the state function 

m 

Sp =—ks inp{x,v,n,m,u:,t) , (5) 

where p{x,v,n,m,uj,t) is the solution of the Fokker- 
Planck equation associated with Q, (0. In our case, Sp 
does not only depend on translational degrees of freedom 
but also includes particle orientation, parametrized by 
n, m and us. The change of particle entropy along a 
trajectory which starts at a point (a:o, vq, Uq, mo, Wq) in 
configuration space at time to and is located at a point 
{x{t),v{t),n{t),m{t),uj{t)) at a later time t is given by 

ASp = kB \np{xo,vo,no,mo,u}o,to) 

- kB \np{x{t), v{t),n{t), m{t),u>{t),t). (6) 

The entropy production in the environment d^env for 
an infinitesimal displacement of the particle is related 
to the heat 5Q released to the thermal bath during this 
displacement, 

doenv — rr ' (^) 
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Following Sekimoto’s stochastic energetics approach [35] , 
the heat is identified with the work done by the particle 
on the thermal environment [40] . This work results from 
the forces the particle exerts on the environment during 
its movement as reaction forces to viscous friction and 
thermal fluctuations, such that 

SQ = -[- 7 U + ovdt 

— [—T]UJ + oujdt ( 8 ) 


for a total particle displacement consisting of a transla¬ 
tional increment da: = vdt and a rotational increment 
ujdt during the time step dt. The symbol o denotes the 
scalar product evaluated according to the Stratonovich 
rule (midpoint regularization) [35]. Note that the heat 
6 Q as well as its two additive contributions from transla¬ 
tion and rotation are scalar quantities and thus invariant 
under rotation of reference frame, a property which we 
exploited in Q by representing the translational part 
in laboratory coordinates and the rotational part in the 
body-fixed frame. The entropy produced by the particle 
in the environment along a stochastic trajectory is the 
integral of over that path with SQ given by Q. Us¬ 
ing (lbI and ( |^ , and the identity [w x (J • a;)] • a; = 0 
we find 


A^env= f I,[f-Vdt' + M-Ujdt' 

Jto ^ 

—mv o dv{t') — (Jo;) o da;(<')] . (9) 

In Appendix [X] we verify that this entropy production 
in the environment can be interpreted as a measure of 
irreversibility [mdH], by relating it to the ratio of prob¬ 
abilities of forward and time-reversed paths for the com¬ 
bined translational and rotational motion. We also verify 
that the total entropy production AS, given by (|4] ) with 
(§, (§ , fulfills the usual fluctuation theorem |38]7 

( 10 ) 

when averaged over the path ensemble with fixed initial 
configuration, denoted by (...). Its average is therefore 
always positive (by Jensen’s inequality) 

(A^)>0, (II) 

in accordance with the second law. 


III. OVERDAMPED LIMIT 

As described in the Introduction, the typical system we 
want to study with our model 0 , (0 is a colloidal el¬ 
lipsoidal particle of micrometer-size (or nanometer-size) 
suspended in water with an inhomogeneous temperature 
distribution around room temperature. It is well known 
that in such systems the effects of the fluid bath in form of 
thermal fluctuations and viscous friction have significant 


influence on the particle dynamics. Indeed, viscous fric¬ 
tion effects are by orders of magnitude larger than inertia 
effects [ini[ii], so that one typically disregards inertia 
completely and adopts the so-called overdamped approxi¬ 
mation by simply putting mass equal to zero, m = 0. In a 
strict mathematical sense, however, this is generally not a 
valid procedure as it changes the order of the differential 
equations of motion 0 , 0 , so that small inertia effects 
(as compared to friction forces) actually correspond to a 
singular perturbation. Therefore, the overdamped limit 
should be performed with some care in a genuine per¬ 
turbative way. Although the overdamped approximation 
and the overdamped limit give identical results in many 
cases, this can not be taken for granted, in particular 
for stochastic differential equations like 0,0 mm, 
involving heterogeneous heat baths, or for functionals 
along stochastic trajectories generated by 0 , 0 , like 
heat [m [45] or the entropy production (0 [^ [26] . 


A. Multiple time-scales 

A standard tool for the rigorous derivation of singular 
limits is the so-called multiple time-scale technique m, 
which exploits that the singular limit is tied to the ap¬ 
pearance of various well-separated dynamical time-scales 
in the system. In our case, the dominance of viscous fric¬ 
tion over inertia effects is related to the existence of two 
distinct time-scales in the equations of motion 0 , 0 . 
The first one is the time over which the velocity degrees 
of freedom relax and reach their stationary distribution. 
It is given by 


m 



( 12 ) 


where 70 is the “typical” friction coefficient of the parti¬ 
cle, e.g., for a spherical particle of radius a it corresponds 
to Stokes friction 70 = diriya {v is the dynamic viscos¬ 
ity of the medium the particle is suspended in). For a 
micrometer-sized particle in water, is of the order of 
microseconds m- The second time-scale is the time af¬ 
ter which one can observe the diffusive motion of the 
particle over a detectable distance, typically about the 
particle size a. 


_ Q^7o 
k^To 


(13) 


with Tq being the average temperature of the bath. Un¬ 
der identical conditions, i.e. for a micrometer-sized par¬ 
ticle in water, is of the order of seconds. Note that 
although we have based our definitions of Ty and on 
translational motion, the time-scales for rotation are ba¬ 
sically the same, because the “rotational mass” (i.e. mo¬ 
ment of inertia) scales as / ~ and the rotational 
viscous friction as ry ~ 70^. 

Based on this separation of time-scales, we calculate 
the overdamped limit of the entropy production 0 [with 
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(§, @] using the multiple time-scale method. It is con¬ 
venient to perform this limit for the generating function 
of the entropy production in the environment @ m- 
As side results, we then find the well-known overdamped 
versions of the equations of motion 0. 0 , as well as 
the relation connecting the probability density p in ([^ 
with its overdamped counterpart. The latter is needed 
to deduce the overdamped limit of the change in particle 
entropy (ro|, and from that the full entropy production 

0. 


B. Generating function 


In order to define the generating function of the 
entropy production in the environment ASenv, we 
first rewrite ® as an integral in dt' only. Observing 
that /cb d[(/aJ) • a;/2fcBT] = /cb d[w^/a;/2fcBT] = 
(lu)) o du)/T — [ci;^/a;/2T^] dT and that also 
/cb d[a;^/a;/2/cBT] = —feBdlniCj^ — (3fcB/2T)dT with 
the Maxwell-Boltzmann distribution for the angular 
velocity 


w 


UJ 


\/det(/) [ Ilj 

'(27rfcBT)3/2 2kBT 


(14a) 


we may write —{Iui)odu}/T in Q as A^b dlni(;j^-f-[(3fcBT— 
Iuj)/2T'^]dT. The term —mv o dt> can be recast in 
a similar way as /cb din -|- [(3fcBT — mu^)/2T^] dT by 
making use of the Maxwell-Boltmann distribution for the 
translational velocity, 


Wy 


/ \ 3/2 

I m\ 


mv^ 

2/cbT 


(14b) 


The dr contributions are a consequence of the inhomo¬ 
geneity of the thermal environment leading to temper¬ 
ature variations along the particle trajectory, which are 
given as dT = (dTjdt + v ■ dT/dx) dt according to (|^, 
where dT/dx = {dT/dxx,dT/dx 2 ,dT/dx^). 

Collecting all these pieces together, we rearrange the 
resulting stochastic integral into three parts [19], AS'reg, 
A^time and A^anomj related to the regular entropy pro¬ 
duction, an entropy production due to time-changes of 
T, and a part of the entropy production related to spa¬ 
tial temperature variations, which we show to yield an 
anomalous contribution in the overdamped limit. The 
final result for ASenv from (|^ thus reads 


AS'env = -A:b lnw„(to)waj(to) T^b lnwy{t)wuj{t) 

+ AS'reg + AStime + ASanom (15) 


with 


ASreg — 


AStime — 


ASa 


'to 


'to 


f ■ V V dk^T M ■ uj 


T 


dx 


T 


dt'. 


, dk^T - Iuj\ dT 
2 T^ + 2T2 ) ' 


a 


5kBT — mv^ 3kBT — Iu:\ dT , , 

V ■ — dt , 


2T2 


2T2 


dx 


(16a) 

(16b) 

(16c) 


By splitting the entropy production into these three 
contributions we can essentially separate the effects of 
time- and spatial variations of temperature, (16b I and 
(16c), from the usual “regular” contribution (16a). We 


point out though that the regular part (16a), contains 


a temperature-gradient term, which is compensat ed b y 
the factor 5/2 (instead of 3/2) in the first term of (16c). 
The appearance of the temperature-gradient is inspired 


I 

by the (a posteriori) observation that the entropy pro¬ 
duction in the overdamped approximation, when inertia 
effects are s impl y disregarded, exactly corresponds to the 
expression (16a). This can be verified by calculating path 
probability ratios, see [TS] and the Supplementary Mate¬ 
rial of [T^ . 

The joint generating function of the three contributions 


(16) to the entropy production is m 


GsiSoSo {x, V, n, m, O), t|a;o, Vo,no,mQ,UJo,to) = (exp(-Si AS'reg - S 2 AStime - Sg ASanom) 

5{x(t) — x)5{v{t) — v)5{n{t) — n)5{m{t) — m)5{oj{t) — u))) , (17) 

where the average is taken over paths with fixed initial conditions a;o, uq, mqi <^o at time to, as before. It can be 
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shown to obey the forward Feynman-Kac formula mm 


dGsiS2S3 

'-^SiS2S3 


f v V dk^T M ■ u}\ /SksT — mv"^ SksT — u)^Iu)\ dT 

* ~Y~ ~ f ' dx T ) I 2T2 + ^ 


+ S3 


/ bk'QT — mv^ Zk^T — uAI lj ^ 

dT' 

2r2 2r2 y 

dx 


G, 


(18) 


Note that for si = 0, S 2 = 0, and S 3 = 0, the generating function represents the probability density p used in (§, 


Goooix,v,n,m,u},t\xo,Vo,no,mo,uJo,to) = {S{xit) - x)S{v{t) - v)S{n{t) - n)5{m(t) - m)5{uj{t) - cu)) 

= p{x,v,n,m,u;,t). (19) 


r 


Accordingly, the operator is the generator of the com¬ 
bined diffusion process for translation and rotation asso¬ 
ciated with Q, ([^; its specific expression is given below 
in Eq. (24) (m dimensionless form). 


We furthermore express / in terms of m and the length- 
scale a, 


C. Dimensionless representation 

Since inertial effects in 0 , (ll) are orders of magni¬ 
tude smaller than friction and other forces, the various 


terms in (18) may be of considerably different magnitude 


as well. For a detailed analysis, we rewrite all quantities 
appearing in (181 by introducing dimensionless represen¬ 


tations of order one, so that the different magnitudes of 
terms show up as dimensionless small (or large) prefac¬ 
tors, which we expect to be related to the ratio of the two 
distinct time-scales Ty and Tx (see (121 and @). Our 
choice is guided by physical intuition and by the charac¬ 
teristics of the systems we intend to model with 0 , 0 . 
Most importantly, due to the separation of time-scales 
Tv ^ Tx, we expect the (translational and rotational) ve¬ 
locity degrees of freedom to equilibrate “instantaneously” 
and become of the order the thermal velocity. In contrast, 
positional degrees of freedom change significantly only on 
the “large” scales a and Tx, so that we measure length 
and time using these units. External forces and torques 
are assumed to be of about the same size as the ther¬ 
mal fluctuating forces. We therefore make the following 
ansatz for relating dimensionfull and dimensionless quan¬ 
tities (denoted by a tilde): 


t = Txt, 


knTo. 

V = \ - V . 


x = ax, 

a 


UJ = 


^bTq 


-u). 


ma^ 

n = h , m = m, 
]\d = k'QTQ^h/E . 


( 20 a) 

( 20 b) 

( 20 c) 

( 20 d) 


I = ma^i, ( 20 e) 

the friction tensors in terms of the “typical” friction co¬ 
efficient 7 o, 


7 = 7o7. 7 = 7oa 77 , 


(20f) 


and the temperature held by the average temperature Tq, 


ksT = /cbTo T ■ 


(20g) 


Plugging the relations (20) into (181 and dehning the 


dimensionless variable Si = k^Si, we obtain the dimen¬ 
sionless form of the forward Eeynman-Kac equation 


- e-^M^ - G,,,,,, = 0, (21) 


with 
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~ d d f d d Q 

^ /i ^ijk^j f H“ ^ ~ 1 (-^ ) i j ^j Q ^ ~\~ {I ^ij^lm^jkl ^k^m 

- 


fiVj _v^l dT^ 
f f\dxi 


5Tvi — VjVjVi 3Tvi — I 


jkUJjiOkVi \ ( df \ 

ySii j 




A7t = - 


S2 


3T - ViVj 

2f'2 


+ 


Sr IijiOj^ujj 

~ 2T^ ' 


dT 

l)i 


( 22 a) 

( 22 b) 

( 22 c) 


where we switched to index notation for convenience with 
summation over repeated indices being understood. In 
( 211 , we defined 


e= J— < 1 
V T^. 


(23 


as 

in 


small parameter, expressing the time-scale separatioi 
m the system. Note that the dimensionless version A"'' o 
the generator of the diffusion process used in (18) i 
given as 

A^ = \£~^d + £~‘^M^ + 


= e-^C^ 


Si— 0 ,S 3—0 


. Si=0,S2=0,S3=0 

-b . 


D. Perturbation expansion 


The dimensionless equation of motion (21) for the gen¬ 


erating function ( |17[ ) constitutes the starting point of our 
analysis, with the goal to derive its overdamped counter¬ 
part in the (singular) limit e = —t 0. To explic¬ 

itly account for the observation that the system exhibits 
dynamics on different time-scales we apply the following 
multi-scale procedure. First, we introduce time variables 
6 and t corresponding to the scales given by Ty and 
(see ( p^ and (13)), and a variable § for the intermediate 
scale TOniifef: 


9 = s '^i, d = £ ^i, T = i. (25) 

We assume that the external time-changes of tempera¬ 
ture, friction coefhcients, forces and torques occur on the 
slow time-scale r only, so that the dimensionless versions 
of ^ read 


T = T{x,t), 

(26a) 

7 = j{x,n,Th,T ), 

(26b) 

II 

(26c) 

/ = f{x,fi,rfi,T ), 

(26d) 

M = M{x, h, m, r). 

(26e) 


Then, we expand powers of e, 

Gsrs,-s^ = + eG(i) + e2G(2) + ... , (27) 

where all the G^*^ are a priori assumed to be functions 
of all variables i, v, n, rfi, uj and of all the three dif¬ 
ferent times 9, d, t defined in (25) (of course they also 
depend on the parameters si, S 2 , S 3 , but we skip the 
corresponding subscripts for notational simplicity). As a 
consequence, the time-derivative in ( 21 ) turns into 


d _2 5 -T 9 d 

dt 99 dd dr' 


(28) 


Inserting (27), (28) into (21) and equating terms of 
equal power in £, we find a hierarchy of coupled equations 
with lowest order (order £“^) 


39 


=0. 


(29) 


From (22b) we see that is an operator in the fast 


degrees of freedom v, u) only, so that their dynamics 
is indeed tied to the fast time 9, as expected. Since 
there are no explicit 9 dependences in (see (26)), and 


since we are not interested in the relaxation processes of 
the fast degrees of freedom on 9 time-scales, we can set 
30^9 jQQ = 0 (for all i) in the following. Then, the first 
three in the hierarchy of equations read 

Ai^G(°) =0, 


(30a) 

7WtG(i) = - £^g(°) , (30b) 

ov 

_yC^tG(2) = _ a7^g(°) + - £^G(i) , (30c) 

3t 3v 

with the first line collecting order £“^ terms, the second 
line order £~^ terms and the third line order £*^ terms. 


The solution to (30a) is 


G^°^ = h, m, I?, t)wv{v\x, T)wa,{99\x, r), (31) 


where Wv{v\x,t) and Wu){9o\x,t) are the dimension¬ 
less counterparts of the Maxwell-Boltzmann distribu¬ 
tions for translational and rotational velocity (at given 


































position and time) defined in (14) [51] , and where 


{x,ri,rh,'d,T) is an unknown function of only the 
slow degrees of freedom according to the general defi¬ 
nition 




{x,n,fh,'d,T) = j diidd) . 


(32) 


To proceed with the higher order equations (30bI, 


(30c) we employ the solvability condition. It states that 


the inhomogeneities on the right-hand sides need to be 
orthogonal to the nullspace of the operator A4 adjoint to 
Aff ( Fredh olm alternative, see, e.g., [^). As we can see 
from (22b), the nullspace of contains the functions 


which are constant in v and u). Therefore, the solvability 
condition for (|30^ reads / didw /dd - = 

0. It is str aightforward to show from the explicit expres¬ 
sion (22a) for that f dbdcDG^G^^^ = 0, so that we 
find g^^^ in (31) to be independent of the intermediate 
time-scale d. 


d-d 


= 0 . 


(33) 


Applying the solvability condition in an analogous way 
to (30c) we obtain 


dg^^^ dg^^'’ 


dr 




(34) 


where we have used J dbdd; = 0. 

We can now state more precisely what it means to de¬ 
rive the overdamped equation of motion for the generat¬ 
ing function This overdamped equation should 

be valid only on time-scales beyond those of 9, after the 
fast degrees of freedom have relaxed and have reached 
their stationary distribution (on 9 time-scales), so that 
they can be integrated out. The overdamped generating 
function is thus given by 


= lim / dbdd) Gg 
£—>■0 J 


( 35 ) 


for which we want to calculate the equation of motion 

dgsiS2S3 


dt 


= lim 

S —^0 


dt 


(36) 


The overdamped limit e = yjTy /r^, —)■ 0 singles out the 
zeroth-order contribution in Gg^gjgg and dGg-^g^Ss/dt. It 


follows from (27), (31) that gsiS 2 S 3 is actually identical to 
g^^\ Moreover, inserting the relations ( [2^ and (28) into 
(36), and using the results we so far obtained from the 
perturbation expansion, namely dG^^^/d9 = 0, as well as 
(Ml and (IMl), we find 


dgsi 


dt 


= / dbdwG^G^^). 


(37) 


In the integral on the right-hand side the d/dvi and 
djduji terms from (22a) evaluate to zero. Our quantity 
of interest therefore becomes 


dgs 


dJv. 


dt 


d 


d 


^ijk [ rx ~ fUg, I t/g) 

dxi V drii drrij ' 


- S3 


fjJii _ MjJcj, 

T f dxi T 

I j j 


- Si 


5fjy.—jy y y. ^TJy. Jf - Q ^ \ dT 


2 t2 


2T2 


with the definition 


J-k(v.Cj)= y dbdd;7r('U,d})G(^^ 


dxi ’ 
(38) 


(39) 


for integrals over (polynomial) function s ttI v, uj) in v and 
t2> multiplying We remark that (38) does not de¬ 

pend explicitly on S 2 , as G*-^^ is independent of S 2 ac¬ 
cording to (30b), with the consequence that in the over¬ 


damped limit the temperature variations with time do 
not show up directly. For evaluating the remaining vari¬ 
ous integrals appearing in (38), it is actually not 

necessary to find the full solution G*^^^ of (30b). In¬ 
stead, they can be evaluated directly from (30b) using 
the known solution ( [M| ) for G*^°\ as shown in Appendix 
IbI The final result reads 
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dg^xi 


dt 


JL 

dXi 




9siS2S3 


^ijk 


9 - d ^ 

-T^nk + 
arij orrii 


iv )jiMi - (fj )jiTepig ( —hq + 

T 


-rUf 


Si 


d h I dT 




dxj \f f dx^ 

/, 1 df\ ff, 1 dT 


d{r^)^ 

dxi 




T Tdxi 


Tdxi 


S 2 S 3 2si 


fi _1 dT\ 
df' 




■s? 


- Si 

irn 

f 


T 


MiMj + Cijkifj )jl I flk 


dMi 

dfii 


ruk 


dMi 

dfhi 


'-MiMjgs^-s 


■ 2,si 


, ~ ^ i ‘^UikUjk 

+ S3(S3-1)^' 


E 


+ -^yhk ) e^jkiv ^)jiMig S,-S2-S3 
UikUjk 


E 


U,kU. 


ik^ jk 


\ df df 


y 3^(fe) ' ^(fe) _|_ 2^(0 ' y-'^(k) 2fj0) /iW j dxidxj 


9siS2S3 , (40) 


IV. OVERDAMPED DYNAMICS 


where we have used on the right-hand side that g^^^ = 

(/gjS 2 S 3 ) and where Uik is defined in Eq. (B51 of Appendix 

0 This expression, together with its interpretation in analogy to ([T^ we obtain the overdamped 

the following Sections constitute the main results of this Umit of the probability density p{x,n,m,t) = 

lim£_j.o / duda;p(a;, u, n, m, a;, t) from gooo- According 
to ( [^ and ( |^ , we can thus read off the overdamped 
equation of motion for p{x,n,m,t) from (40) by setting 

Si = 0 , S 2 = 0 , S 3 = 0 : 


dp 

dt 


•^ovcrP ^ 5 


(41) 


with 


At = _ 

ror 


d_ 

dxi 


(7 hfj - (7 ) 


dx 


k^T 


^ijk 


dn. 


rik 


drrii 


rUk 


iv - (77 ^)jik^Tepiq 


dn„ 


drrij, 


(42) 


where we switched back to dimensionful quantities. Here 
and in the following, we use a sub- or superscript “over” 
to indicate the overdamped limit and to differentiate 
overdamped quantities from the original ones. The 
Langevin equations, which are equivalent to the Fokker- 


Planck equation (41), (42), read 


i = 7 V 


7 ^ dk^T k^T dj ^ 

2 dx 2 dx 

+ \/2A:bT7"^^^ o 1(0 : 


(43a) 


n = 
rh = 


g-^M+y2yifg-^/^C{t) 

g-^M + y2kBTg-^^^Cit) 


X n , 
X m , 


(43b) 

(43c) 
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where ^ [t] and C (t) are Gaussian white noise sources (like 
in 0 , §, but not related to those), and where all the 
products involving noise terms, even the cross products, 
are to be interpreted in the Stratonovich sense. These 
Langevin equations are well-known as a model for the 
overdamped translation and rotation of ellipsoidal parti¬ 
cles [T]. In particular the “splitting” of the translational 
diffusion tensor D = into k^T and 7 “^ between 


We now compare the expressions (45) and (46), ob¬ 


tained from applying stochastic thermodynamics prin¬ 
ciples “naively” to the overdamped equations of motion 
(43), with the full result (40) from the systematic pertur¬ 


bation analysis. We first observe that the overdamped 
limit e —>■ 0 of the probability density p appearing in 
the change of particle entropy ^ is given by eval¬ 
uated at Si = 0 , S 2 = 0 , S 3 = 0 , which we obtain from 


the two spatial derivatives in the diffusion term of (42) (31) to be p{x,n,m,t)wv{v\x,t)waj{i^\x,t) (dimension- 


has been obtained before , as well as the specific “spuri¬ 
ous drift terms” in (43a) resulting from that “splitting” 
Note that the explicit form of these “spurious 
drift terms” and the corresponding interpretation of the 
noise products are a natural outcome of the systematic 
perturbation analysis. They can not be obtained from 
the overdamped approximation, naively disregarding in¬ 
ertia effects in 0 , 0 . The overdamped approximation 
yields correct equations only for homogeneous thermal 
environments when T, 7 and rj are constant. 


ful units). Therefore, lim^-^o A^p yields the terms listed 
in (45), plus additional terms involving logWyiv^j which 


cancel precisely with the first line in (15). We conclude 


that, in the overdamped limit, the change of the parti¬ 
cle entropy can indeed be identified with (45), while the 


entropy production in the environment originates from 
the overdamped counterparts AS'°™(jj of 


the entropy terms in the second line of (15). These over¬ 


damped entropy contributions are encoded in our main 
result (40) by the definition 


V. OVERDAMPED ENTROPY PRODUCTION 
AND ANOMALOUS ENTROPY 


Based on the overdamped Langevin equations (43) we 


can follow the standard reasoning of stochastic ther¬ 
modynamics and define the entropy production along 
stochastic trajectories as 


^gover ^ ^gover 


AS° 


(44) 


It is composed of contributions from the entropy change 
of the particle AS'°™"' and from the entropy change in the 
environment A5'™®''. The entropy change of the particle 
is given as [38] 

A^'over ^ lnp{xo,no,mo,to) 

— k^\n.p{x{t),n{t),m{t),t), (45) 

for a trajectory which starts at at a point {xQ,nQ,mQ) 
at time to and is located at a point {x{t),n{t),m{t)) at 
a later time t. The entropy change in the environment 
can in principle be defined from the heat exchanged with 
the environment. However, due to the variations of tem¬ 
perature with position such an identification is subtle for 
the translational degrees of freedom |T9l [49| and thus the 
definition of entropy production in the environment is 
better based on path probability ratios [T5]. It reads 


f 1 

A cover / 

‘^*^env / fji 

Jto ^ 


t 


(yy-^M) ■Mdt' + o dW(t') 


(46) 


where dW is the increment of the Wiener process corre¬ 
sponding to the Gaussian white noise C (^) in (43b), (43c). 
In that way, one arrives at an entropy production which 
is expressed as a sequential functional on overdamped 
trajectories. 


ffsiS2S3 n,m,t\xo,no,mo,to) 

= (exp(-siA5— - S2A5,7“ - ssASZl^J 

S{x{t) — x)5{n{t) — n)5{m{t) — m)) , (47) 


analogous to ( [17| ). To access and analyze the specific form 
of an individual contribution we set the s variables asso¬ 
ciated with the other contributions to zero and compare 
the remaining terms in (40) with the general formulas 
from Appendix [0 In that way, we find the following: 

The entropy production given in (46 1 arises precisely 


from the overdamped limit of the regular part (16a), rep¬ 


resented by the Si terms in (40). There is no contribution 
from the entropy production (16b) due time-changes of 


temperature, because (40) does not explicitly depend on 


S 2 (as already pointed out earlier), such that is 

bound to vanish identically. However, the S 3 terms yield 
additional contributions to the overdamped entropy pro¬ 
duction, which are not included in ( |46[ ). From their spe¬ 
cific functional form we infer (cf. Appendix [0 that these 
contributions can not even be expressed as a sequential 
functional over overdamped trajectories. Their origin is 
the entropy production AS'anom from (16c) (see also 0 ) 


Although it is not possible to explicitly write this 
“anomalous entropy production” [Tn| as an integral along 
paths of the overdamped dynamics, we can still derive 
a number of interesting and useful results on its aver¬ 
age behavior from (40), similarly to the reasoning in the 


Supplementary Material of m- Setting Si = 0, S 2 = 0, 
S 3 = 1 (i.e. Si = 0 , S 2 = 0 , S 3 = l/fee) we find that gooi = 
{exp{-ASZoni/kB)S{x{t) - x)S{n{t) - n)5{m{t) - m)) 
obeys the same forward equation as 5000 = n, m, t), 
so that it has the solution 


(exp(-A5°rm/fcB) 

S{x{t) — x)6{n{t) — n)S(m{t) — m)) 

= p{x,n,m,t) ■ (48) 
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Integrating over the spatial coordinates a;, n, m, we find 
the fluctuation relation 

(exp(-A5romAB)) = l. (49) 

It follows immediately (by Jensen’s inequality) that 


The explicit form for the average rate of anomalous en¬ 
tropy production can be obtained from ([40|) by observing 


that = -f^Jdxd{n,m)dgoosJds 3 \ 


dt \—"^anom/ dt J \ ^ - -o i qq — ii 

(see also the derivation of Eq. (CIO) in Appendix [C| . 
It reads (in dimensionful quantities) 







d 

dt 


(A5aTo)n) 


kBi 


E 


2T /A ^(fe) _|_ 2-y(0 ^ ^(fc) _|_ 27^(0/!^ 


E 


dT 




dT 


\ 2T 


E7E, + E 


kB / da; d(n, m) 


j^fdT\ 


7 -f 27*'*^I 

2 


E 


V I 


2 r7(') 

id) jm 


-IN 


2T \dxJ 




■+i:(^ 


-1 




dxj 

dT dT 
dxi dxj 
-1 


dT 

dx 


(51a) 

(51b) 

(51c) 


These expressions for the average rate of anomalous en¬ 
tropy production are another central result of the present 
paper. We give three different, but equivalent forms. In 
the first line, the entropy production is written in the 
coordinate frame fixed to the particle. Accordingly, the 
quantity Uik^^ is the temperature gradient along that 
principal axis of the particle, for which the friction coef¬ 
ficients are 7 *^^^ and 7 *^^^ (and moment of inertia is 
see Eq. (B5) in Appendix [b] where also the Utk are de¬ 
fined. The second line represents the anomalous entropy 
production in the laboratory frame of reference (I de¬ 
notes the identity matrix). In the third line, we switch 
back to vector notation and we express the average (...) 
over particle trajectories explicitly as an integral over the 
probability density p = p{x,n,m,t); note that the inte¬ 
gral over d(n, m) is not performed independently over n 
and m but rather represents an integral over the space 
of particle orientations, parametrized by n, m. 


The result ( [sT] ) generalizes the discovery of the anoma¬ 
lous entropy in |19j in essentially two respects. Eirst, it 
covers non-trivial particle shapes, quantifying deviations 
from a perfectly spherical bead by the principal values 
7 ^), r]d) and id) jm. Note that the ratio id)/m does 
not depend on the particle mass, as id) is proportional 
to m; it is thus a purely geometrical factor reflecting 
the particle’s shape. Second, it takes into account the 
rotational degrees of freedom of the Brownian particle, 
and reveals that rotational motion adds to the entropy 
production. Partial results about the effects of rotation 
on the entropic anomaly have also been derived in |26j . 
The rotational entropy production originates from the 
terms involving the rotational friction coefficients rjd) ^ 
since these terms vanish when we “freeze” the particle 
rotation using the limit rjd) —n 00 to be left with the 
translational motion only. We can therefore identify the 
translational and rotational contributions to the anoma¬ 


lous entropy production as 


K».n. = ^(?7-‘ + E(^ + 2r“’l) 


_ kBP 

^rot — 2^ 

I 



Id) I'm j 



, (52a) 


(52b) 


It has been argued in m that the anomalous entropy 
is generated by the particle permanently transporting 
heat between adjacent regions at different temperatures 
in the inhomogeneous thermal environment on the fast 
time-scale Ty and associated length-scale (feeTowj^/^/yo 
without performing any visible displacement on the long 
(overdamped) time-scale We can therefore interpret 
Ktrans and Krot as State-dependent, anisotropic heat con¬ 
ductivities quantifying this process behind the anoma¬ 
lous entropy production |26] . This interpretation also 
explains why the rotational contribution Krot depends on 
the translational friction coefficients 7 . The conducted 
heat is “stored” in the rotational degrees of freedom, but 
is transported from one temperature region to another by 
translation. Without translational motion, the particle’s 
rotation can not produce “anomalous entropy”. 

We finally remark that the rotational motion con¬ 
tributes to the anomalous entropy production even if the 
particle is perfectly spherical. For homogeneous spher¬ 
ical beads with radius a, translational and rotational 
friction tensors as well as the moment of inertia ten¬ 
sor are proportional to the unit tensor, i.e. we have m 

^(1) = -y(2) = ^(3) _ jy(l) = 7y(2) = ^{ 3 ) — 

and id) = id) = id) = 2md} so that the conductivi- 
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ties become isotropic, 
S/cbP 

^trans — /is — 


^rot — 


6T7(b 

SkBP 


5kBP ■ 
iQ-nTva 
1 


I = 




2T ^(1) _|_ 2 r;(i)/^ 927rTz/a 


(53a) 
I.(53b) 


For the translational part we just recover the result from 
m, while the particle rotation gives rise to an additional 
contribution not described in [TS]. The frictional “coef¬ 
ficients” 7 ^^^ and representing translation and 

rotation, respectively, quantify physically related effects 
and are of similar magnitude. Hence, the rotational con¬ 
tribution to the “anomalous entropy” is well comparable 
to the translational part, and is actually only by about a 
factor four smaller (see also the discussion in [ 2 H] 1 . 


geometry (and medium viscosity), although the result 
would be quite lengthy and cumbersome and is explicit 
only up to the quadratures in (56 1 . As it turns out. 


however, these integrals can be performed analytically in 
the case of spheroids, i.e. ellipsoids with two equal semi¬ 
axes, ai = 02 - One distinguishes between oblate and 
prolate spheroids. 


Oi = 02 > 03 oblate spheroid, 
oi = 02 < 03 prolate spheroid. 


A. Flat oblate spheroid 


(58) 

(59) 


For oblate particles, as defined in (58 1 , we obtain from 


(56) 


VI. PROLATE AND OBLATE SPHEROIDS 


For an ellipsoidal particle 


2 2 2 

Jb-\ JjO JUr> 

^ + ^ + ^ = 1 
of of of 


(54) 


with semi-axis lengths oi, 02 , 03 , the translational fric¬ 
tion coefficients 7 ^*^ have been calculated by Oberbeck 
[S3], the rotational ones 77 ^®^ by Edwardes [M] and Jef¬ 


fery I 
read 


They are summarized, for instance, in |34] and 


Ol = Q!2 = 


2 

0^3 3 

o? 


a.i 

arccos ( 

03 

ai 

1 V 

“1/ 

1 - 4 

(i-S 

^3/2 

ai 

arccos ( 

03 

03 

V 

“1/ 


X = 


l_ll| / ai\3/2 

(S) 


arccos 


ai 


(i-s) 


1/2 ■ 


(60a) 


(60b) 


(60c) 


with 


7 ^*^ = 167rj2- 


X- 


afai 


(for z = 1,2, 3), 


7 ( 1 ) = 


77 ( 2 ) = 


77 ( 3 ) = 


X = 


IGmy 

of -f of 

3 

afa 2 + o-sO-s 

IdTTl/ 

of -I- of 

3 

a§a 3 - 1 - of«! 

IdTTl/ 

of -I- of 

3 

of Oi -I- of 02 


dA 

A(A)’ 

dA 


Qlj — 


0 (of + 


A(A) = [(of + A)(of + A)(of 


(55a) 

(55b) 

(55c) 

(55d) 


(56a) 

(for i = 1,2,3), (56b) 

+ A)]1/2, (56c) 


and with 1 / being the viscosity of the medium. The mo¬ 
ment of inertia tensor for the ellipsoid (54) is 


For further analysis we focus on the limiting case of a flat 
oblate spheroid, assuming 


Oi 


(61) 


Neglecting second and higher order terms in <5, the trans¬ 
lational and rotational friction coefficients then read 


^(1) = ^(2) ^ 1^1 + y(5 

7 ( 3 ) = 16z/ai 

32 


77(3) = |z7a3(l + A^ 


4 
— ( 
TT 


(62a) 

(62b) 

(62c) 

(62d) 


The moment of inertia tensor has only zeroth and second 
order terms in S, so that to first order we find simply 


/(i) 0 0 

0 0 I = 

0 0 /(3) 


/ 2,2 

/ ° 2+°3 

5 

0 

0 


0 

5 

0 


0 

0 

a^+al 

5 


(57) 

It is clear that from these expressions we can write 


down the entropic anomaly (511 in terms of the particle’s 


1 


/(I 


0 0 


^00 


— 0 /(2) 0 = 0 ^ 0 
' 0 0 /(3) I \ 0 0 ^ 

5 


(63) 


From these expressions we can easily calculate the heat 
conductivities (52) associated with the anomalous en- 
































13 


tropy production to first order in S, 

(64a) 
(64b) 
(64c) 
(64d) 

where the = 15 2,3) are the translational heat 

conductivities along the principal axes of the particle de¬ 
fined via (Ktrans)ij = ^^i/ct(,fcKtrins 5 and likewise for the 
rotational conductivities K(ot- We find that both, trans¬ 
lational and rotational contributions to the anomalous 
entropy production are comparable in magnitude, with 
the translational one being about four times larger than 
the rotational one. 


^trans 


^(3) 

'^trans 


.( 2 ) 


k-BP 

32Th'a^ 


kBP 

32Tvai 

13 


19 

~8 


67 

T^' 


7 


di) 


( 2 ) 


- 


,(3) 


kBP 

32Tva^ 


kBP 
?>2Tvai 
147 


64 
497r 
23 

44 ~ 21787r' 
120 


2201 


1697r 


B. Thin prolate spheroid 


For the prolate particles from (59), we can write (56) 


as 


Otl — 0,2 — 

ai 


arccosh ( ^ 


(S) 


2 

^3 3 


1 - 


-1 


3/2 


(1-4) 

arccosh 


X = 


('-4) 

arccosh 


3/2 


03 


(- 4 ) 


1/2 ■ 


Considering thin prolate spheroids with 

— =5 < 1, 

03 


we obtain to lowest order in <5 


7 ( 1 ) = 7 ( 2 ) = 


(3) _ 




21n2-hl-21n(5’ 


7 ' 1 = 


^( 1 ) = ^( 2 ) = 


8711203 


21 n 2 - l-21n(5 ’ 
167rj2a| 

3(2Tn2“l”“2lnJ) 


d3) _ l^_,,„3x2 


rj' ' = 




(65a) 

(65b) 

(65c) 

( 66 ) 

(67a) 

(67b) 

(67c) 

(67d) 


These friction coefficients asymptotically vanish in the 
limit (5 —> 0 (with a logarithmic approach to 0 , except for 
77 *-^^), because then the prolate particle more and more 


resembles a one-dimensional rod-like object which expe¬ 
riences practically no friction when moving through the 
fluid. It is easy to verify though that even for 5 —>■ 0 
inertia effects remain negligible compared to viscous fric¬ 
tion forces, since the particle mass, being proportional to 
the particle volume, decreases much faster with (5 —>■ 0 
than the friction coefficients (67). In other words, the 
condition t^/tx ^ 1 [see (23)]as a prerequisite for the 


overdamped limit is fulfilled for arbitrarily small 5. The 
moment of inertia tensor //m as well has entries which 
vanish asymptotically as <5 —>■ 0 so that we here have to 
keep the second-order terms. 


/ Id) 

0 

0 \ 

/ 

5 

0 0 \ 

1 ° 

/( 2 ) 

0 I = 


0 

ald+sd 0 

V 0 

0 

/(3) J 

[ 

0 

0 ^52 y 


(68) 

Calculating the heat conductivities (Ktrans)z/ = 


fdjfcCjfeKarins and (Avrot)z7 = from (1^ USffig 


(fc) 


the expansions (67), ( 68 ), we find that the leading order 
terms diverge with ln(T/5). Explicitly, the results read 


di) 


( 2 ) 


^trans ^trans 


kBP 


32Tz/a3 V 37r 5 


II, 1 Illn2-b4 
In ^ 


37r 


k-Bp 7^28 1 2(70 In 2-19) 


= 

ra na 


(69a) 

“ 32Ti2a3 ~5 25^r ' ’ 

( 1 ) ^ ( 2 ) ^ kBP fl2_ 1 3(2601n2 + 99) ^ 

rot «rot ^ 2 Tl 2 a 3 V1377 6 84577 

(69c) 


( 3 ) kBP ^ 24 1 , 3(401n2 + 3) 

^rot — ~~~ m -I- 


S2Ti>a3 \2377 6 


11577 


(69d) 


Comparing with the corresponding results (641 for the 


flat oblate spheroid, we see that the anomalous entropy 
production rate typically is larger for the thin prolate 
particle, due to the logarithmic divergence of the con¬ 
ductivities. 


VII. SLIGHTLY DEFORMED SPHERE 


Another interesting case to consider is an ellipsoid 
with almost identical semi-axes, i.e. a particle slightly 
deformed from perfect spherical shape. We fix ai and set 


02 = ai(l-I-(52), a3 = ai{l + S3), 


(70) 


assuming 


^2 < 1, (53 < 1 • 


(71) 


Plugging (70) into (56), we can now perform the integra¬ 
tions by expanding the integrands in ^ 2 : < 53 , and calculate 
the friction coefficients from (55) to a desired order in 62 , 
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63 - The first-order results read 
7 ^^^ = 67ri/ai 
= dirvai 


2 2 

1 -I- -62 + -S 3 
5 5 


( 2 ) 


7 


7 ( 3 ) = dirvai 


. Ir 2^ 

1 -I- -62 + -S 3 
5 5 


■,(1) 


= Sttz/cV 


. 2^ 1 


= Sttz/Oi (1 -I- |(52 -I- 1^3 
?y(3) ^ SiTval (\ + ^S2 + 1^3 


(72a) 

(72b) 

(72c) 

(72d) 

(72e) 

(72f) 


Obviously, the well-known isotropic Stokes friction coef¬ 
ficients for translation and rotation iszi of a perfectly 
spherical particle are recovered in the limit S 2 —> 0, 
Js —> 0. 

For the anomalous heat conductivities Ktrans: ^rot [see 
(52)] we further need the moment of inertia tensor //m, 
which is easily obtained by inserting (70) into (57). We 
finally find, again to first order in S 2 and i53, 


di) 


72) 


'^trans 


71) 


^(2) _ 
'Irot ~ 


73 ) 


5fcBP 

SdirTvai 

^kBP 

SdirTvai 

5kBp 

SdirTvai 

3fcBP 

927rTi/ai 

3fcBP 

927rT^'ai 

SkBP 


75 75 

75 ^ 75 ^ 

75 ^ 75 ^ 


, 118, 118, 

345 ^ 345 ^ 

, 109, 118, 

345 345 ^ 


1 - 


118 , 109 , 

-O 9 — -Oq 

345 345 ^ 


(73a) 

(73b) 

(73c) 

(73d) 

(73e) 

(73f) 


92TTTvai 

For (^2 = 0, ^3 = 0, we recover the results for a spherical 
bead, already calculated in (53). Note that our observa¬ 
tion from the spherical case, that translational and rota¬ 
tional contributions to “anomalous” entropy production 
are well comparable (about a factor four difference), also 
applies for near-spherical particles. Deviations from per¬ 
fect spherical shape lead to similar corrections for both, 
translational and rotational “anomalous” entropy. 


VIII. CONCLUSIONS 

In the present paper we analyze the thermodynamic 
properties of a single Brownian particle of non-spherical 
shape, for which rotational degrees of freedom play a 
non-negligible role. The main aim is to understand how 
the stochastic thermodynamics is affected by “coarse- 
graining” the level of description of the particle motion 


when performing the overdamped limit to “integrate out” 
the fast velocity degrees of freedom. This question is 
of particular interest in case that the surrounding heat 
bath is heterogeneous, a situation which is known to be 
non-trivial already for the particle’s equations of motion 
[4TM48j . A central quantity for such an analysis is the 
trajectory-wise entropy production of the particle de¬ 
fined according to stochastic thermodynamics [16j . For 
translational Brownian motion, it has been discovered 
in |19j that the overdamped limit of this entropy produc¬ 
tion generates an “anomalous” contribution, which is not 
captured by the statistics of the overdamped trajectories. 

Here, we analyze in detail the effects of a non-spherical 
particle shape and of the Brownian rotation of such par¬ 
ticles. Starting from the standard entropy production 
of stochastic thermodynamics (extended to include rota¬ 
tional degrees of freedom, see Sec. II) on the level of the 
full-fledged description of the particle dynamics, we per¬ 
form the overdamped limit using singular perturbation 
theory (Sec. III). As our main result we find that the ro¬ 
tational Brownian motion not only yields a “standard” 
contribution to entropy production which is consistent 
with the overdamped approximation (where one simply 
disregards velocity degrees of freedom), but in addition 
also generates an “anomalous” entropy which can not be 
expressed as a functional along overdamped trajectories 
[see Eqs. (51) and (52)]. This “anomalous” contribu¬ 
tion to entropy production from the particle’s rotation 
is comparable in magnitude to the “anomalous” entropy 
generated by translational motion. 

We remark that our starting equations ([^, ([^ to 
model the particle’s Brownian motion do not contain any 
hydrodynamic coupling between translational and rota¬ 
tional degrees of freedom. Such couplings would be rele¬ 
vant, for instance, for particles with a helical shape. For 
that reason our analysis is restricted to the class of el¬ 
lipsoidal (and other rod- and disk-like) Brownian parti¬ 
cles. We expect that the more general case of hydro- 
dynamic couplings between translation and rotation will 
also induce couplings between these degrees of freedom 
in the entropy production and thus lead to an additional 
“anomalous” contribution. The details, however, remain 
to be revealed in future work. 

Potential applications of the present findings include 
the influence of rotational Brownian motion and the as¬ 
sociated “anomalous” entropy production in inhomoge¬ 
neous thermal environments on optimal time-dependent 
protocols, realized by external forces to optimize a spe¬ 
cific quantity of interest during a finite-time process 
on the efficiency of microscopic stochastic heat engines 
[28l [56] , and even on its universal fluctuations discovered 
recently in m- 
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Appendix A: Entropy production 


In Sec. |TT]of the main text, we derive the entropy pro¬ 
duction in the environment ([^ from the heat exchanged 
between particle and thermal environment. In this Ap¬ 
pendix, we summarize the relevant results which show 
that ([^ can also be obtained from the ratio between the 
probability for observing a certain particle trajectory and 
the probability for observing its time-reversed counter¬ 
part [TB] . 

We consider a set of stochastic equations of the form 


Mi = , (Ala) 

dpi = {-TijPj + Ui) dt + Pij o dWj , (Alb) 

where qi are space-like coordinates, transforming as qi —>■ 
qi if time is reversed, and pi are velocity-like coordinates 
with Pi —>■ —pi under time-reversal. The dynamics of the 
space-like coordinates is linear in the pi with a tensor /i^ 
which may depend on qi only (not on pi). The determin¬ 
istic part of the dynamics for the velocity-like coordinates 
Pi has a dissipative contribution —TijPj and a part Ui col¬ 
lecting external forces, which is assumed to transform as 
Ui —>■ Ui under time-reversal. The tensor Pij defines the 
strength of the thermal noise, where the stochastic noise 
itself is represented by the increments of mutually inde¬ 
pendent Wiener processes dW^; the products between /3y 
and dWj are to be interpreted in Stratonovich sense. The 
diffusion tensor Dij resulting from these noise terms is 
given by Dij = PikPjk, where we assume that its inverse 
{D~^)ij exists. Note that the structure of the stochastic 


differential equations (A1) covers the Langevin equations 
Q and (|^ for translational and rotational Brownian mo¬ 
tion. 

It can be shown (see, e.g., [15] or the Supplementary 
Material of IISI) that the probability P[q{t),p{t)\qQ,pQ] 
for a specific trajectory {q{t),p{t)) starting at {qQ^Po) at 
time to and ending at the point {qi,Pi) at a later time 
ti is related to the probability P[q{t),p{t)\qQ,pQ] for the 
time-reversed trajectory [i.e. with (qgjPo) = {Qij~Pi) 
and (qi,Pi) = {qo, -Po)] according to 


^[g(0,P(0|go:Po] 

P[Q{t),P{t)\Qo,Po] 


= exp 



2(T1 )ijPikPk [Uj 


Pj) 



(A2) 


We observe that for subsets of coordinates with noise 
sources which are statistically independent, the inverse 
diffusion tensor D~^ takes block-diagonal structure, so 
that the contributions of these subsets in the exponent 
are additive, even though their deterministic forces may 


depend on the whole set of coordinates. As the Gaus¬ 
sian noise sources in Q and ([^ are independent, we can 
therefore focus on translation and rotation separately to 
determine their contributions to the exponent in (A21. 

For the translational motion Q we identify q = x and 
p = V and find by comparison with (|A1|) 


Pij — dk 


r — 


— fi 1 
m 

Pd — ■ 


Id 
m ’ 


(A3a) 

(A3b) 

(A3c) 


The exponent in (A2) thus reads 
2{D~^)..TikPk{uj - Pj) - d^dpi = (uj/j - 
TTiViVi)/(k^T) and corresponds exactly to the transla¬ 
tional contribution in up to a factor 1 /fee. 

Likewise, the rotational motion (§ matches the equa¬ 
tions (Al I for q = (n, m), p = lj and 


Pd 

r 

Ui 

Dij 


f ^ijk^k 0 1 

\ 0 eijkrrik ) ’ 

(A4a) 

^ik‘nkj 7 

(A4b) 

(/ (/ )ij^jkl^k^lm^m 7 

(A4c) 

2 kBT{r^),k{I-%Vki- 

(A4d) 


We can then calculate its contribution to the exponent 
to read 2 TikPk (uj - pj) - dui/dpi = (w^M* - 

lijUJiCoj)!(ksT), which is equivalent to the rotational part 
in ([^, again up to a factor \/k^. 

Summarizing, we therefore arrive at the central result 
of this Appendix, 


^ _ p-AS„,/fcB 

P~ 


(A5) 


relating the entropy production in the environment 
AS'env as defined in 0 to the ratio of forward and back¬ 
ward path probabilities. Finally, from here it is straight¬ 
forward to verify |38j that for the total entropy produc¬ 
tion from Q the integral fluctuation relation 


(A6) 


is fulfilled, as already stated in (10). The average (...) 
is taken over all trajectories starting from a fixed initial 
condition. 
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Appendix B: Evaluation of the integrals 




We here describe how to evaluate the integrals 

Jvi = J (Bla) 

=y’dhdcDwiG^i), (Bib) 

JvjVjVi = J dvdu}VjVjV^G^^\ (Blc) 

/ dudw/,fcWjWfe0,G(i), (Bid) 


appearing in (38), without calculating G^^^ expli citly . 
The basic idea is to multiply the order e~^ equation (30b) 


by a polynomial n('i;, lj), integrate over i3, cD, and rewrite 
the left-hand side in terms of the operator Ad adjoint to 
(see ( |22b[ )). In that way we obtain equations of the 

form 

J dudcD (Adn)G(i) = - j dudw n(£l'G(°)), (B 2 ) 

where we remembered that 9G ^°^/dd = 0 according to 
(33). For suitable choice of 11, the operation Adll in the 


left-hand integral may then reproduce one of the desired 
polynomials from (Bl| (or a linear combination of sev¬ 
eral such polynomials), while the right-hand side can be 
calcu lated straightforwardly using the exp licit expression 
(22a I for and the known solution ( |3li for G*-®^. 


For instance, taking 11 = {1^ we get Adll = such 

that after evaluation of the corres pon ding right-hand side 
integral —f dvdd) Vi(£^G^^^) in (B2| we find 


Jv.- 1)(7 ')b 

(B3) 

Similarly, with the choice 11 = we obtain 


df\ 


J 


= -iv )ikT€jki ( 


d 


d 

drhj 


m g 


AO) 


- (si - l)(f? 


(B4) 


The two remaining integrals (Blc) and (Bid) are 
more cumbersome to calculate because of the third- 
order polynomials involved. For Jy-v-j,- it turns out 
that we have to choose n = AijuiVjVkVi, where 
the tensor Aijki has to be determined such that in 
= -AijklijjmVmVkVl + JkmVmVjVl + ^Ira H" 

2 f AijkiijjkVi + jjiVk + %iVj) the sum over polynomials 
of third degree in ii-components reduces to VjVjVi. To 
construct the explicit form of Aijki, it is convenient to 
diagonalize 7 in Ad by 

{U^Wh = , (B5) 

where t/ is a symmetric tensor and 7 ^*^ are the eigenval¬ 
ues of 7 . It is then straightforward to identify 




UimUjmUkn Ui n 


j(m) _|_ 2^(") 


(B 6 ) 


Note that the sum here is over m and n, and that 
Aijki obeys the symmetries Aijki = AjM = Aijik- 
With this expression for Aij^i we can now calcu¬ 
late the two integrals / dhdcD {MAijkiVjVkVi)GA') and 
— / dvdiv AijkiVjVkVi{£^G^°^) to find 


JvjVjVi — 5TJv. -|- (s 3 — 1)T 


‘^UikUjk 


+E 


3^1'^) 

U.kUjk \ df 


^ j{k) _|_ 27(0 J dxj 


(0) 


(B7) 


The sum over I is specified explicitly, since this index ap¬ 
pears only once; apart from that summation over dou¬ 
ble indices is still understood. To arrive at the sim¬ 


ple form (|B7|) we made use of (|B3|) and of the re- 
Ukll = 3 ^UA~~ ~ l^l 

i^dA^kij d~ 2Aijki^^ki d~ {2 Aiik i -(- Aikud^ki — bdii, which 


lations 2 Aiiki + Ami = 


can be proven by using (B5), (B 6 ). 


Finally, the calculation of <7/.^-.-^-. from (Bid) pro¬ 
ceeds completely analogously. The proper choice for 11 
turns out to be BijkiVjCjkOJi, with 


Bijkl — 


U^mUJ^VknVlJ^^'> 

/IM 


(B 8 ) 


to guarantee that MBijkiVjOJk^i = Ijk'^i ^kV i + 
(terms linear in h-components). In determining (B 8 ) we 


have assumed that I and fj are diagonalized simultane¬ 
ously by V [58j , 

{V^IV),j = . (B9) 

Then, evaluating ^ dvdCb {MBijkiVjCjk(^i)G^A and 
— / dvdu) BijkiVjCjkCbi^DG'A^) we obtain 




(S3-I)f^ 


U,kU. 


ik^ jk 


dT 


j{k) _|_ 2fjd)/id) dxj 


,( 0 ) 


(BIO) 


Again, we have simplified ( |B10[ ) by using (B3), and by 
observing that Bijki{i~^)ki = - 'Ei 
‘2^Bijkl{^I ')km{d )^nf?r) 


' Biklm{,d — dSij. 


Appendix C: Forward equation for the generating 
function of sequential functionals 

We consider the general Langevin-equation in Ito form 

dq^ = Uidt + ■ dWj , (Cl) 

where the dot here denotes the Ito product. The set 
of coordinates Qi typically comprises velocities and posi¬ 
tions for translational motion, but, more generally, may 
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also contain angular velocities and corresponding coordi¬ 
nates representing the particle orientation. Note that 
we therefore adopt a slightly different notation as in 
Appendix where we explicitly distinguished between 
space- and velocity-like coordinates. 

The deterministic velocities are ui, the dWi are in¬ 
crements of independent Wiener processes, and the j3ij 
define the symmetric diffusion tensor Dij via Dij = 
PikPjk- In the general case, Ui and jdij are functions of 
qi and t. The Langevin-equation can also be written in 
Stratonovich interpretation. 


where the trajectories start at < 7 ( 0 ) = go nt time to nnd 
end at q{t) = g at a later time t. The functions /i, gi, 
fi depe nd on coordinates g^ and time t. The products 
in ( |C4| ) are understood in the Ito sense. The equivalent 
form of the functional with Stratonovich products is 

J{q,t\qf),to) = f [hdt'+giodqi{t') + fiOdWi{t')] . 
Jto 

- - 

Like in ( |C2[ ) we label the the functions h, gi, fi of the 
Stratonovich form by an overbar. They are related to the 
Ito functions via 


dgi = Uidt + Pij o dWj , (C2) 

where the relation between Ui and Ui is 

Ui = Ui + ^ ■ (C3) 

2 dqk 


h — h — Dij + - — pji 

2 dqj 2 dqj 

, (C 6 a) 

9i ~ 9i ■> 

(C 6 b) 

fi = fi- 

(C 6 c) 

The generating function of J{q,t\qQ,tQ) 

is defined as 


The general form of a functional along trajectories gen¬ 
erated by (Cl I reads 


J{q,t\qo,to) 



' +gt- dqift') + fi ■ dWi{t’ 


)], 

(C4) 


Gs(q,t|qo,to) = (e - g)) , (C7) 

with the average (...) being taken over all trajectories 
starting at go, to- It can be shown (see, for instance, 
Chapter 6.4 in [T] or the Supplementary Material of [12] 
for the explicit derivation in special cases) that the for¬ 
ward equation for the generating function Gs reads 


dGs 

dt 


A^Gs 


—s {h + utgi) Gs + 


2 


{Dijgtgj + fifi) Gs + s 


A 

dqi 


{Dijgj -(- Pijfj) Gs 


(C 8 a) 


r 


with the usual generator of the forward diffusion process state g, we obtain 


= -p- (ui - p-Dij 

OQi \ OQj 


(C 8 b) 


dt ' ' 

= - / f sti -b smg^ - pD,jg,gj - pftfi ) e 


-sj 


(C9) 


All functions Dij, Pij, Ui, h, gi, fi in (C 8 ) are evaluated 
at the final state g, t. 


Using the definition (C7) and integrating over this final 


Finally, deriving with respect to s and setting s = 0 
afterwards, we find an equation for the average of the 
functional ff over all trajectories starting at go, Iq, 


^ (J) = {h + u,gi) . 
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